Griffiths-McCoy singularities in random quantum spin chains: Exact results 
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We consider random quantum (tight-binding, XX and Ising) spin chains in the off-critical region 
and study their Griffiths-McCoy singularities. These are obtained from the density of states of 
the low-energy excitations, which is calculated exactly by the Dyson-Schmidt method. In large 
finite systems the low-energy excitations are shown to follow the statistics of extremes and their 
distribution is given by the Frechet form. Relation between the Dyson-Schmidt technique and the 
strong disorder renormalization group method is also discussed. 



I. INTRODUCTION 

Quenched disorder has a profound effect on the low- 
energy, low-temperature and long wavelength properties 
of quantum systems. The interplay between quantum 
fluctuations, correlations and disorder fluctuations gen- 
erally results in strong singularities in the thermody- 
namical quantities and in the (dynamical) correlation 
functions 1 -^. This type of effect takes place even out- 
side the quantum critical region, e.g. in the quantum 
paramagnetic phase at zero temperature, T = 0, where 
spatial correlations are short ranged^ 3 -. The origin of 
this phenomenon, as pointed out by Griffiths 4 , is due to 
rare regions, in which strong bonds are accumulated by 
extreme fluctuations, so that the system in these regions 
is locally in the thermodynamically unstable ferromag- 
netic phase. As a consequence the excitation energy, E, 
in the rare regions is very small, the relaxation process is 
very slow and the associated relaxation time, r ~ E , 
is divergent in the thermodynamic limit. If we consider 
a finite part of a sample with linear size, £, the charac- 
teristic time-scale of the slowest relaxation process stays 
also finite and asymptotically given by: 

r-r. (i) 

Here z = z(S) is the dynamical exponent, which is gener- 
ally a continuous function of the quantum control param- 
eter, S, which measures the distance from the quantum 
critical point. 

According to scaling theory 5 -! 6 - the distribution of the 
low-energy excitations, n(E,£), depends on the scaling 
combination, E£ z , and for a small but fixed E it is pro- 
portional to the volume, £ d , since the probability of find- 
ing a rare region goes linearly with the volume. From 
this the asymptotic behavior of the distribution function 
in the thermodynamic limit reads as: 

n(E) ~ E d ' z - X . (2) 

Thermodynamical quantities which are obtained through 
integration of the density of states are also singular. For 
example the low-temperature behavior of the average lin- 
ear susceptibility, x{T), and that of the specific heat, 



c v (T), is expected to scale as^ 3 - 

X(T) ~T- 1+d / z , c v {T)~T d ' z , (3) 

whereas the small-field, H, dependence of the zero- 
temperature magnetization is given by: 

m(H) ~ H d/z . (4) 

One can see from Eq.Q that the susceptibility is diver- 
gent at zero temperature for z(S) > d, which was noticed 
first by McCoyZ. in an exact calculation of the random 
transverse- field Ising chain (RTFIC). 

Detailed results about Griffiths-McCoy singulari- 
ties are obtained for one-dimensional (Id) systems 
partially by numerical investigations (free-fermionic 
techniques'^, density matrix renormalization 
metho d 10 ' 11 , quantum Monte Carlo (MC) simulations 12 ) 
and by analytical calculation s 11 ' 13 ' 14 ' 15 based on 
the use of a strong disorder renormalization group 
(SDRG) method 2 . For higher dimensional systems 
Griffiths-McCoy singularities are studied numerically, 
either by quantum MC simulations 16 or by numerical 
implementation of the SDRG method^. 

Analytical and conjecturedly exact results about 
Griffiths-McCoy singularities are scarce and these are 
practically restricted to the RTFIC. Analytical solution 
of the SDRG equations is first obtained in the vicin- 
ity of the quantum critical pointi 3 -, i.e. in the weakly 
disordered and in the weakly ordered Griffiths phases, 
where the dynamical exponent is shown to diverge as 
z(8) ~ 1/|#|- The solution is then extended to the com- 
plete Griffiths phas e 11 ' 15 and the calculated value of z(5) 
is shown to agree with that obtained through a mapping 
to a random walk problem in a random environment* 1 ^ 

In this paper we use a direct and simple method to 
calculate exact values of the Griffiths-McCoy singulari- 
ties in a class of random quantum spin chains. These 
models include the random tight-binding chain, the ran- 
dom antiferromagnetic XX-chain as well as the RTFIC. 
The low-energy excitations for each models have the same 
form: they are obtained from the eigenvalue problem of 
a symmetric tridiagonal matrix, M. see in Eq.0, with 
random (positive) entries. In the off-critical region of the 
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spin chains there is an even-odd asymmetry: the matrix- 
elements of M. are taken from different distributions at 
even and odd bonds. We calculate the density of states, 
n(E), in the center of the band by the Dyson-Schmidt 
technique^ using the random walk idea by Eggarter and 
Riedinger— . In Ref.[2(| n(E) is calculated in the contin- 
uum approximation for an even-odd symmetric M. , which 
corresponds to the critical point of the random quantum 
spin chains. In the present paper M. has a general even- 
odd asymmetric form, which corresponds to a strongly 
disordered quantum Griffiths phase and for which the 
continuum approximation is no longer valid. Having the 
exact behavior ofn(E) at hand we than calculate the sin- 
gularities of the thermodynamic quantities (specific heat, 
susceptibility, magnetization). 

The structure of the paper is the following. Random 
quantum chain models studied in this paper are presented 
in SecHIl In Sec lIIII the density of states of the low energy 
excitations is calculated by the Dyson-Schmidt technique 
and the relation of this technique with the SDRG method 
is discussed. Thermodynamic singularities are calculated 
in Scc|lV] and the results are discussed in SecfVl 



II. RANDOM QUANTUM CHAINS 

A. Random tight-binding model 

The first model we consider is a one-dimensional tight- 
binding model with off-diagonal disorder 21 being defined 
by the Hamiltonian: 



w=5Xii)<<+i| + |t + i><i|) 



(5) 



with random hopping matrix-elements, i,-. The hopping 
matrix-elements are generally taken from different distri- 
butions at even (t e ) and odd (t ) sites, so that a quantum 
control-parameter is defined as: 

[lnt ]av - [lni e ] 
var(lnt e ) + var(lnt ) ' 

where [. . . ] av stands for averaging over quenched disorder 
and var(x) stands for the variance of x. For 5 > (5 < 0) 
the model is asymmetric and the particles arc preferen- 
tially at odd (even) bonds. The symmetric model with 
(5 = corresponds to a quantum critical point. 

In the basis, the Hamiltonian is represented by a 
tridiagonal matrix 



M 



[0 h 
h 

ta tg 



(7) 



and we arc interested in its eigenvalue problem: 

Ma = Ea (8) 

and the corresponding density of states, n(E), at the 
center of the band. 



B. Random antiferromagnetic XX-chain 

The second model is the random antiferromagnetic 
XX-chain defined by the Hamiltonian: 



Hxx = ^ Ji(SiSi +1 + SfSV +1 ) 



(9) 



in terms of the spin-1/2 operators, S^ y , at cite i. Here 
the Ji exchange couplings are random variables which 
have different distributions at even (J e ) and odd (J D ) 
sites. Using the Jordan- Wigner transformation: = 



± iS? and c[ — af exp iri 



-l + - 



and 



exp 



a 3 a 3 



a i this Hamiltonian is expressed in 



terms of fermion creation (c^) and annihilation (cj) op- 
erators as—: 



7~Lxx 



i c i c i+l 



+ h.C.) 



(10) 



The low-energy states of the model contain one fermion, 
which can be written in the form = ^2iCticf\0), 
where |0) denotes the fermionic vacuum. Energies in this 
one fermion subspace are obtained by the solution of the 
eigenvalue problem of M. in Eq.([7]) with the correspon- 
dence: 



(11) 



Then the quantum control parameter of the model is just 
given by S in Eq. ^ . In the asymmetric model with <5 > 
(5 < 0) there is enforced dimerization and the system is in 
the random dimer phased with preference of odd (even) 
bonds. On the other hand at the quantum critical point 
with (5 = the system is in the so called random singlet 
phase— 



C. Random transverse-field Ising chain 

Our third and final model is the RTFIC, which is a 
prototypical model of random quantum systems having 
an order-disorder transition— This system is defined by 
the Hamiltonian: 

i i—l 

in terms of the Pauli-matrices, a*' z , at site i and the Ai 
couplings and the hi transverse fields are random num- 
bers. 

As for the XX-chain Hi is expressed in terms of 
fermion operators 2 ^: 



1 



c c% 

, 2 



Hi = -J2 h > 

i 

- \Y, x ^°t -^t+i+ci+i) (1 3 ) 
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which is than diagonalized through a canonical transfor- 
mation. Now the low-energy excitations contain one free 
fermion, the possible energy of which is given by the pos- 
itive eigenvalues of the following symmetric matri^iS^: 



T 



/0 hi 
hi Ai 
Ai 

V 



(14) 



This is equivalent to M in Eq.f7|) with the correspon- 
dences: 



t 2 i-l 



to 



A; 



(15) 



Using this relation together with Eq.© the control pa- 
rameter of the RTFIC is given by the difference in the av- 
erage log- fields and the average log-couplings. For 5 > 
(5 < 0) the system is in the paramagnetic (ferromag- 
netic) phase, and 6 — represents the quantum critical 
point. 

We can thus conclude that the low-energy properties of 
all the three models are related to the eigenvalue prob- 
lem of M. in Eq.([7]). In the next section we calculate 
the density of states of matrix M. around E — by the 
Dyson-Schmidt method. 



III. DENSITY OF STATES AT THE CENTER 
OF THE SPECTRUM 

Here in the first two subsections we recapitulate the ba- 
sic ingredients of the Dyson-Schmidt method and present 
the solution in the continuum approximation. Our 
findings, which are obtained in the strongly disordered 
regimes are presented in the last two subsections. 



A. The random walk method 

In order to calculate the density of states of M. 
we introduce a new vector, A, with the components 
A, = on-iti-i/ ai, which satisfy the equations : Aj+x = 
tf/(E — Ai). The basic ingredient of the Dyson- 
Schmidt method^ is the node counting theorem of one- 
dimensional Hamiltonians, which states that the inte- 
grated density of states, N(E) = J_ n(E')dE', is given 
by the fraction of positive terms in the sequence of A,. 
At the center of the band, E = 0, the components 
of A have alternating signs, thus here the "sign vari- 
ables" Si = sign[Ai(— l) 4 ] have a fully ordered state, 
. . . iiiiiimf ... an d jv(0) = 1/2. For nonzero E 
the iterated equations for Ai are the following^: 



A2i — [ , 2% ) A 2 ;_ 2 



/2i-2 — 



^2i-2 , 

1 - £/A 24 _ 2 



1 + (£A 2l _ 2 - E*)/tl_ 2 



(16) 



which lead to different iteration behaviors for small posi- 
tive E for various limiting values of A 2 i . These are sum- 
marized as: 



A 2l+ i/A 2l < 0, 

hi = i, 

A 2l+2 /A 2l < 1, 



if 
if 
if 



A 2l < E (17a) 
E > A 2l > P/E (17b) 
A 2l « P/E (17c) 



where t denotes the typical (average) value of the matrix- 
element. According to (|17bj) we can identify an interval, 
[E,P/E], in which the "signs" stays ordered, say =|. 
There is a finite upper boundary value at A max = P /E, 
where the iterated sequence is reflected, but Si stays f 
(see (|17cj) . whereas as the sequence arrives at the lower 
boundary value, A min — E, the "spins" change sign (see 
(|17ap and the iteration process starts again, however in a 
new domain with Si =|. Consequently for a small E > 
the sign variables have a fragmented domain structure 
•■■4444TTTTT4444--- an d therefore the fraction of posi- 
tive terms in the sequence of Aj is somewhat larger than 
1/2, due to extra positive terms appearing at the domain 
walls. If the typical (average) size of a domain is denoted 
by £, then the density of states is asymptotically given 
by: 



N(E) - N(0) = 



1 

2l 



(18) 



We can thus summarize that to obtain the density of 
states at the center of the spectrum it is enough to fol- 
low the evaluation of the sequence, Aj, within one typi- 
cal domain and calculate its size, t. Within this domain 
we formally put / 2 i- 2 — 1 in Eq. (|16p and set i) a re- 
flecting boundary at A max and ii) an absorbing bound- 
ary at A m i n . If we introduce the logarithmic variable, 
lnA 2 j = U2i, we obtain a random walk (directed poly- 
mer) problem: 



u 2 i = 2(lni 2i _i - lni 2l _ 2 ) + u 2 i~2 



(19) 



with reflecting (u = u max ) and absorbing (u = u m i n ) 
boundary conditions. In this language the walker (poly- 
mer) starts at uo = u max and its mean first-passage time 
(length) at the position u m m is just £, thus u m i n = Uc ■ 



B. Analysis in terms of the diffusion equation 

In order to set the length-scales in the random walk 
problem we use a continuum approximation in which 
Eq. (|19p is transformed into a diffusion equation: 



dP(u,£) _ D d 2 P(u,£) 



dP(u,£) 
du 



(20) 



Here P(u, £) is the probability distribution of the walk, 
D = 2[var(lnt e ) + var(lnt )] is the diffusion coefficient 
and v — 2([ln£ ] av — [ln£ e ] av ) is the drift velocity. The 
typical size of the transverse fluctuations of the walk is 
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given by: u = D/v = S^ 1 , whereas the average distance 
between two reflections, £, follows from the relation u ~ 
\fD£, thus we obtain for the correlation length: 

£,~D- l 5- 2 (21) 

which agrees with the result of SDRG calculations^. The 
continuum approximation and thus the use of the diffu- 
sion equation is justified if the correlation length is much 
larger than the lattice spacing. This condition is satis- 
fied if we are either at the critical point, 5 = 0, or in the 
weakly disordered Griffiths phase with \S\ <C 1. 



the band there is a power-law singularity of the density 
of states: 

AT(£)-jV(Q)~(£) ' . (26) 

which is equivalent to the form in Eq.([2]). This result for 
the random antiferromagnetic XX-chain has been pre- 
sented in2S. 



C. Analysis in the strongly disordered Griffiths 
phase 



1. Critical point 

At the critical point both the correlation length, £, and 
the typical size of transverse fluctuations, u, are divergent 
and they are related with the length scale, £, as: £ ~ £ 

and u ~ \J Di. Absorption of the walker in this case is 
due to typical fluctuations, when u grows to the order 



of the width of the strip: u 
ln(P/E 2 ). From this follows: 

1 



Au 



so that 



-HP/e 2 ) 2 



N(E) - N(0) ~ D [\n(t 2 /E 2 )] ' 



(22) 



(23) 



This is the classical result derived by Eggarter and 
Riedinger— . 



2. Weakly disordered Griffiths phase 

In the weakly disordered Griffiths phase with 1 -C S > 
the walker is drifted towards the reflecting boundary 
and both the correlation length and the typical size of the 
transverse fluctuations are finite, but much larger than 
the lattice spacing, thus the continuum approximation 
is valid. In this case u is much smaller than the width 
of the strip, its absorption takes place with a very small 
probability: p(Au) cx exp(— ^p), thus it is a rare region 
effect and due to extreme fluctuations. Before having 
such a large fluctuation the walker is reflected several 
times and the typical number of independent excursions 
is given by £/£, the value of which follows from extreme- 
value statistics^: p(Au)£/£ = 0(1). From this we have 



£exp(-^ln(P/£ 2 



l/z 



with 



1 - 2v - 2 
z D [varQnt, 



[hlt ]av - [lute 



var(hi t )] 



= 26 



(24) 



(25) 



Here z is just the dynamical exponent defined in Eq.([T]). 
In the weakly disordered Griffiths phase at the center of 



In the strongly disordered Griffiths phase the correla- 
tion length is in the order of the lattice spacings and the 
continuum approximation is not valid. In this case we 
use discrete variables and denote the (nonlogarithmic) 
position of the walker at the j-th step of the fc-th inde- 
pendent excursion, which starts at r(k), as A^. Thus 
r(k)/k — £ for large k and the normalized position is 
given by: 



(fe) _ 
= 



,(fc) 



j+r(k) 

n 

3'=l+r(fe) 



to 



L 2j'-2 

The condition of absorption is formulated as: 

A 



(27) 



(k) 

min min p\J 

k l<j<Ar(k) J 



E 2 



^ > (28) 



where Ar(fc) = r(k + 1) — r(k), which can be replaced by 

Ar = oo. Keeping in mind that is typically much 
larger than its minimum value we can estimate the order 
of magnitude of the minimum as 



mm mm p 2 

k Kj/<oo 



(fe) 



,(fe) 



(29) 



Here ?/ fc ' is a Kesten variable^ for any k, the distri- 
bution function of which for large arguments displays a 
singularity: 



p(y) 

where the exponent, 
the equation: 



-(i+aO 



(30) 



/i, is given by the positive root of 



= f 



(31) 



(For a pedagogical introduction to the theory of Kesten 
variables see Appendix C of Ref.0.) In this way the 
typical number of excursions, £/£, follows from extreme- 
value statistics^: £/(, f°° p(y)dy — 1, and we obtain: 



£. 



2fi 



(32) 
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Comparing with Eq. (j24]) we see that the dynamical ex- 
ponent in the strongly disordered Griffiths phase is given 
by: 

- = 2/i , (33) 

z 

which in the limit 5 <C 1 gives back the result obtained 
in the weakly disordered Griffiths phas o 11 ! 15 in Eq. (f2"5} . 
Then with the correspondence in Eq. (j3"3"|) the density of 
states at the center of the band is given in Eq. ([26|) . 

D. Relation with the strong disorder 
renormalization group method 



TABLE I: Analogous quantities in the random walk (RW) 
and in the SDRG methods 



method 


independent unit 


length scale 


energy scale 


RW 


excursion 


size of the excursion 


(k) 


SDRG 


cluster 


size of the cluster 


t'o(k) 



to an excursion (between two reflections) of the walk of 
size Ar(fc) ~ £ and the minimal value of p£) for this 
excursion is just the renormalized value of t' {k). The 
analogous quantities in the two approaches are collected 
in Table [J 



The density of states in the center of the spectrum of 
M. can be analyzed by the SDRG method^, too, and here 
we outline this procedure. The first step in this study is 
to arrange the matrix-elements, U, in descending order 
and use the largest one, Q = maxi{^}, to set the energy 
scale in the system. Let us denote the largest term by tj , 
which connects sites j and j + 1 and eliminate the two 
equations in the eigenvalue problem which contain tj. 
In second-order perturbational method, which is correct 
up to O ({tj-i/tj) 2 ) and O ((tj+i/tj) 2 ) we have for the 
effective matrix-element, t', between the remaining sites, 
j — 1 and j + 2: 



t' 



tj-itj+i 



(34) 



This new term has a length, m! — rrij-i +rrij +mj+\ = 3, 
where the original matrix-elements have unit lengths. 

In the following steps we repeat the decimation trans- 
formation, during which the energy scale is reduced, the 
lengths are increased and the distribution functions of 
the matrix-elements, R e (t e ,fl) and R (t ,Q), approach 
their fixed-point form. This type of RG equations have 
been analytically solved both at the critical point^ and 
in the Griffiths phas o n i 15 . Here we summarize the known 
results for the Griffiths phase with S > 0. 

In the starting steps of the RG both t e and t Q terms are 
decimated, but the transformation in later steps become 
asymmetric. As the typical lengths are growing beyond 
ml ~ £ almost exclusively the t' terms are decimated and 
the t' e terms become very small, such that at the fixed 
point, O — > VL* = 0, we have t' e /t' — ► 0. As a consequence 
the energy of the low-energy excitations is simply E — t' Q . 
At the fixed point the distribution of t Q is given by^i^: 



2fi 



(35) 



where fi is defined in Eq. (|31[) . This is just equivalent to 
the distribution of the excitation energies in Eq.@, with 
the dynamical exponent defined in Eq. ([33|) . 

Now to make a correspondence with the random walk 
method the starting RG steps which lead to an effec- 
tive t' e (k)(^> t' {k)) of length m'(fc) ~ £ are equivalent 



IV. THERMODYNAMIC SINGULARITIES 

Here we consider the random tight-binding model with 
half filling, as well as the random antiferromagnetic XX- 
chain and the RTFIC and note that all these models are 
expressed in terms of free fermions. The common form 
of the Hamiltonians is given by: 



H F = Y,E q {r)+r, q - 1/2) 



(36) 



where E q denotes the g-th eigenvalue of M. and 77+ 
(r) q ) are fermion creation (annihilation) operators. The 
ground-state energy per site of this system is given by: 



2L ^ 9 2 



and the free energy per site: 



n(E)EAE 



(37) 



E„ 



2T 



- T I In 2 



Er, 



n(E) In 



cosh ( 

2T 



dE U38) 



where L is the length of the chain. From the free energy 
we obtain the internal energy: 



S(T) 



1 



E„ 



n{E)Et&n\i[^)<lE >:]•)• 



and the specific heat: 



E„ 



»W ( w 



cosh 



^-) NO! 



Now using the form of the density of states at the center 
of the band we obtain for the low temperature behavior: 



/oo 
e 1 /^ 1 coslT 2 e de. , 
-OO 



(41) 



in agreement with the scaling result in Eq. ([3]) . Note that 
the prefactor in Eq. (|4"l"T) , A, is proportional to £ _1 z~ 1 , 
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TABLE II: Relation between the length-scales in different 
regimes of the quantum control parameter. 



critical point 


5 = 




weakly dis. Griffiths 


\5\ « 1 


e> £ > a 


strongly dis. Griffiths 


|5| = 0(1) 


I > i ~ a 



which means that in the weakly disordered Griffiths 
phase we have: A ~ S 3 , 6 <C 1, in agreement with the 
SDRG result 13 . 

Next we consider the random antiferromagnetic XX- 
chain for which in the Hamiltonian in Eq. @ we introduce 
a homogeneous ordering field: HJ2i^ z - This term with 
fermionic variables assumes the form: H/2Y2A c t c i — 
1/2), thus the eigenvalue matrix M. contains also diag- 
onal elements: A4i t i — -ff/4, Vi, and the eigenvalues are 
shifted by E — > E + H/4. The magnetization is obtained 
through differentiation: 

m(H,T) =--^~J n(£)tanli ^— J dE (42) 

where we have used the fact that the spectrum of M. 
in Eq.([7j) is symmetric to E = 0. At zero temperature 
m(H, 0) is singular for small H: 

m(H, 0) - N{H/A) - N{-H/A) ~ H 1/z (43) 

as in Eq.(j4]). Evaluating the integral in Eq.(j42|) for small 
H and T, however with H/T = 0(1) we obtain for the 
low-temperature susceptibility: 

X (T) ~ T 1 ^- 1 (44) 

which corresponds to the scaling result in Eq.([3]). 

V. DISCUSSION 

In this paper we have studied Griffiths-McCoy singu- 
larities in random quantum (tight-binding, XX and Ising) 
spin chains, which can be represented in terms of free 
fcrmions. The main step of our investigation is the cal- 
culation of the density of states of the low energy exci- 
tations, which excitations are eigenvalues of a symmetric 
tridiagonal matrix with random entries, however with an 
odd-even asymmetry This latter problem is solved ex- 
actly by the Dyson-Schmidt techniqu e 19 ' 20 for any value 
of the quantum control parameter, 5. Previous studies 
of this problem are restricted to the quantum critical 
ponrb 2 ^, 8 — 0, and to the weakly disordered Griffiths 
phased, <5<1. 

As we described in Sec lIIIBl in this problem there are 
three length scales: the mean-first passage length, i, the 



correlation length, £, and the lattice spacing, a. In the 
different regimes of the quantum control parameter their 
relative magnitudes are summarized in Table HTl 

In a finite system there is still another length scale 
given by the size of the system, L, and the mean- 
first passage length can not exceed this value: £ ~ L. 
Consequently the lowest excitation energy is limited to 
Ei ~ L~ z . In this case one is interested in the distri- 
bution of the scaling combination, E\L Z , which in the 
random walk method in Sec lIII CI is obtained from the 
statistics of extremes. Here we recall that E\ is just 
the minimum value of a set of L/t; independent random 
numbers, each having the same parent distribution in a 
power-law form, see Eq. (|3"0)) . Consequently the distribu- 
tion of ei = aE\L z in the large L limit follows the Frechet 
distribution^!: 

PiM^e^expt-e^) , (45) 

where a is a nonuniversal constant which depends on the 
amplitude of the tail in Eq. (f3"01) . Here one can go on 
and consider the second eigenvalue, E2, or more generally 
the q-th smallest eigenvalue, E q . These are all obtained 
from the theory of extreme value statistics of indepen- 
dent and identically distributed (i.i.d.) random numbers 
and their distribution is given by the generalized Frechet 
distribution, see ir>21. In this way we have shown that the 
distribution of the lowest energy levels of these strongly 
correlated physical systems are described in a form which 
holds for i.i.d. random numbers. This scenario, which 
is shown here exactly for the specific models is expected 
to hold generally for all such random quantum systems, 
even in higher dimensions, for which the low-energy be- 
havior is controlled by a so called strong disorder fixed 
point in the SDRG framework^. 

The dynamical exponent, z, which is calculated exactly 
in this paper is found a continuous function of the control 
parameter, 5. Using the SDRG approach the same re- 
sult is obtaine d 11 ' 15 , thus our present study gives further 
credit to the conjecture that the SDRG method provides 
asymptotically exact results even far outside the critical 
point, as far as dynamical quantities are considered. This 
latter statement is expected to hold for all systems with 
a strong disorder fixed point. 
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